How the small hyperfine splitting of P-wave mesons evades large loop corrections 
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The recent discoveries of the bottomonia states /i6(lP) and hb{2V) confirm the quark model 
prediction, already verified in the charmonia sector, that the hyperfine splitting of P-wave mesons 
is very small. The striking agreement is somewhat surprising because the nonrelativistic result, for 
which the splitting is zero, may be modified due to large mass shifts from coupling to open flavour 
meson pairs. This paper is based on the observation that in most models hyperfine splitting remains 
small despite what are in many cases large mass shifts. This effect is shown to be a generic feature 
of models in which the coupling is driven by the creation of a light quark pair with spin-one. 
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Quark potential models share a common spin- 
dependent structure, which in perturbation theory 
yields an expression for the mass Mgi^j of a state of 
spin 5, orbital angular momentum L and total angu- 
lar momentum J, 



MsLj = M 
AsL.j = As' 



A 



SLJ, 



(1) 
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in terms of common expectations values M, As, At 
and Ao of the spin- independent, spin-spin, tensor and 
spin-orbit terms. For states with L ^ 0, certain linear 
combinations of masses are independent of three of 
the four expectation values, and the most interesting 
of these is the hyperfine splitting. For the P-wave 
family it is expressed: 
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(Af3p„ + 3Af3p^ 4- 5M3pJ - A/ip^ = A,,. (2) 



In the nonrelativistic limit A^ = 0, and the experi- 
mental charmonia masses are in excellent agreement 
with this prediction [1] : 



M 



Xc(lP) 



/lc(lP) 



-0.02 ±0.19 ±0.13 MeV. 



Recently the BaBar collaboration discovered the 
hbilP) in the decay T{3S) /ibTr+Tr" [2], while Belle 
observed the same state, and discovered its radial ex- 
citation /ib(2P), in the process e~^e~ — >■ /ifeTr+Tr" [3]. 
The corresponding splittings are also very small: 



^Xb(lP) ~ ^46(1P) 
^Xb(lP) - ^45(1?) 



M 



Xb(2P) 



2±4± 1 MeV(BaBar), 
1.62 ± 1.52 McV(Belle), 
MeV(Belle). 



Lattice QCD calculations likewise exhibit very small 
hyperfine splittings [4-6]. 

That the hyperfine splitting is so small is a triumph 
of the nonrelativistic quark model, but it is also some- 
thing of a surprise. Relativistic effects, as well as ad- 
mixtures of different orbital or radial components in 
any of the wavefunctions, will cause deviations from 
the nonrelativistic result. Moreover, it may be antici- 
pated that mass shifts due to coupling to open flavour 









AM3p, 


AMip^ 


Ind. 


BS(lP,cc) 


459 


496 


521 


504 


-L8 


K (lP,cc) 


198 


215 


228 


219 


-1.3 


LMC (lP,cc) 


35 


38 


63 


52 


-2.9 


YLCD (lP,cc) 


131 


152 


175 


162 


-0.4 


OT (lP,cc) 


173 


180 


185 


182 


0.0 


OT (IP, 66) 


43 


44 


45 


44 


-0.4 


OT (2P,66) 


55 


56 


58 


57 


0.0 


LD (1P,66) 


80.777 


84.823 


87.388 


85.785 


-0.013 


LD (2P, bb) 


73.578 


77.608 


80.146 


78.522 


-0.048 



TABLE I. The magnitudes of the mass shifts computed 
in various models. The final column "Ind." shows the 
induced hyperfine splitting due to loop effects. 



meson pairs will lead to further discrepancy. The aim 
of this paper is to show that the latter effect is small. 

The effect of "unquenching" is to shift the masses of 
physical states downwards with respect to their bare 
"quenched" masses. Owing to their different spin and 
total angular momentum quantum numbers, states of 
given multiplet (such as xo-, Xi: X2 and h) have dif- 
ferent couplings, and the resulting mass shifts cause 
deviations from the quenched mass formula (1). 

Barnes and Swanson [7] derive a loop theorem in 
the limit in which mesons sharing the same flavour, 
orbital and radial quantum numbers have equal mass. 
Remarkably, all such states are shifted by the same 
amount, so that the effect of unquenching can be ab- 
sorbed into a redefinition of model parameters. In 
this limit the equal mass shifts cancel exactly in equa- 
tion (2) and the nonrelativistic result of zero hyperfine 
splitting is maintained. 

The same authors consider the more realistic case 
of mass shifts due to loops of D, D*, Dg and D* with 
physical masses, using also physical masses for the 
various charmonia. In principle there will be further 
corrections due to the different masses in loops with 
excited mesons, such as Dq, Di, D'j^ and D2, but it is 
reasonable to expect that these are smaller since they 
are suppressed by an energy denominator. The mass 
shifts AMsLj calculated in this way are shown in the 
first row (BS) of Table I. 
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Although the relative shift between any two of the 
states is small (« 50 MeV) compared to the overall 
mass shifts (« 500 MeV), it is still very large com- 
pared to the scale of the experimental hyperfine split- 
tings. It is thus striking to note that the correction to 
equation (2) due to unquenching, 



-^(AM3p„ 



- 3AM3p^ + 5 AAf3p^ ) + AAfip^ , (3) 



is just -1.8 McV. 

Could this be accidental? For comparison Table I 
presents the mass shifts and induced hyperfine split- 
tings of charmonia and bottomonia IP and 2P states 
in various other unquenched quark models: those of 
Kalashnikova (K) [8], Li, Meng and Chao (LMC) [9], 
Yang, Li, Chen and Deng (YLCD) [10], Ono and 
Tornqvist (OT) [11], and Liu and Ding (LD) [12]. Al- 
though the models differ markedly in the magnitude of 
the shifts, they share the common feature that the in- 
duced hyperfine splitting is in each case significantly 
smaller than relative mass shifts among the states, 
which are themselves significantly smaller than the 
overall mass shifts. The same, however, cannot be 
said of states which lie above open flavour threshold, 
and all of the calculations that follow apply only to 
subthreshold states. 

The effect is particularly noteworthy given the dif- 
ferent assumptions underlying the various models. 
While the shifts in the approach of BS are derived in 
second order perturbation theory, those of K, LMC, 
YLCD, OT and LD are obtained by solving the 
coupled-channel equations. While BS, K and LD use 
harmonic oscillator wavefunctions (either of universal 
size or with different sizes for charmed and charmed- 
strange mesons), LMC, YLCD and OT use wavefunc- 
tions obtained by solving a coulomb plus linear poten- 
tial model. In BS, YLCD and OT, the pair creation 
strength is flavour-independent, while for K, LMC, 
and LD the creation of strange quarks is suppressed 
with respect to that of light quarks. 

One feature common to all models is that the cou- 
pling is driven by the creation of light quark pair in 
spin triplet. The quark spin and spatial degrees of 
freedom factorise so that the amplitude for the cou- 
pling can be expressed as a linear combination of spa- 
tial matrix elements, which are the overlaps of the 
meson spatial degrees of freedom, weighted by angu- 
lar momentum recoupling factors. For a state with 
S, L and J quantum numbers coupling to a pair of 
S-wave mesons with spins Si, S2, there is a single spa- 
tial matrix element Ai for each partial wave / [13]. 
The corresponding recoupling coefficients Cglf can 
be derived from the general expression of ref. [14] and 
for the case of P-wave mesons are shown in Table II. 

The magnitude of the mass shift due to a given 
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^SLJ 
^IID 




1 




1/4 
3/4 


3/20 
9/20 
2/5 



1/2 
1/2 


2-/S1S2 ^SLJ 


1/2 


1/4 


-1/4 






TABLE II. Angular momentum recoupling factors. 



channel is 



dp 



(4) 



while the probability (in perturbation theory) that the 
physical state is in the corresponding two-meson chan- 
nel is 



^SLJ 



'-■SLJ 



dp 



p^\Mp)\ 



(5) 



The energy denominators are written in nonrelativis- 
tic form and are functions of the reduced mass and 
binding energy. 



^SLJ 



ms2 



SLJ, 



(6) 
(7) 



for loop mesons with mass rus^ and • If one uses 
the physical mass for Mslj then (4) is a coupled- 
channel equation; if instead one uses the bare mass 
for Mslj then (4) is the second order perturbation to 
the energy shift. 

The total mass shift and continuum probability are 
the sums over the corresponding quantities for the dif- 
ferent spin channels and partial waves: 



AM. 



SLJ 



E ^^^s\ 

si S2I 



Sr>l 

SLJ ' 



P' 



SLJ 



EpSiS2l 
^SLJ 

S1S2I 



(8) 



It is reasonable to assume that mesons which share 
the same orbital and radial quantum numbers but dif- 
fer in quark spin and total angular momenta have 
identical radial wavefunctions. Thus, for example, D 
and D* have the same radial wavefunctions, as do xo, 
Xi, Xi s-iid ^- In that case the spatial matrix element 
Ai is independent of the spin and total angular mo- 
menta, and the dependence on these quantum num- 
bers only enters in the energy denominator and the 



si S2I 
SLJ 



coefficient C 

In the limit in which the loop theorem of Barnes and 
Swanson [7] is derived, all channels are characterised 
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by a common reduced mass /.i and binding energy e, 
and everything can be written in terms of common 
integrals, 



AM' 



dp 



(9) 

Since the integrals are common to all channels the 
summations (8) involve only the recoupling factor, 
and from Table II, 



SLJ 



1, 



(10) 



which leads to the loop theorem. 

To account for departures from the equal mass limit, 
consider a situation in which the bare masses of the 
quarkonia states are given by the perturbative formula 
(1), while those in the loop are given by corresponding 
formula 



m + (5(ii)si, 



(11) 



and are characterised by centre-of-mass m and split- 
ting 6. 

To find relations among the mass shifts, one can 
set up a power series expansion in which everything is 
expressed in terms of spin-averaged quantities, which 
correspond to setting all spin splittings to zero (/z = 
to/ 2 and e = 2m — M). For a given channel the re- 
duced mass and binding energy are related to their 
spin-averaged counterparts by 



MsiS2 



2m ^ ^ ' \TO^ / 

^SLJ-^ + H(^-2)s. + {^k)s2)~^SLJ. (13) 

Defining X^^^} by 

M.,..4L} = M<l + ^5Lj)> (14) 
one can express the mass shift as a power series 



AM 



S1S2I _ QS1S2I t^SiS2 1 



SLJ 



SLJ 



X E^-XSLJY' 
n=0 



dp 



p'\Mp)\ 



(l-fpV2/ie)" 



+1 ■ 



(15) 



The series converges due to the smallness of Xg^^j and 
the monotonic decrease of the integrals with n. The 
resuhs of BS [7], K [8] and LD [12], with whose decay 
amplitudes the integrals above can be computed ana- 
lytically, are very accurately reproduced keeping only 
the first three terms in the expansion. Keeping instead 
only the first two terms (which serves as a reasonable 
approximation for cc and an excellent approximation 



for bb) the mass shift can be expressed in terms of the 
spin-averaged values AAf' and of equation (9), 



A ]\ fSlS2l ^ /^SiS2l ^SiS2 

^^^'^SLJ ~ ^SLJ ~~~ 



(AM'-X^i^j6P'). (16) 



Ignoring terms suppressed by higher powers in the 
small parameter (5/ to, the mass shift is 

(17) 

where 



y-lS\S2l /^S\S2l 

^SLJ — ^ 



5LJ ((H>«i + (H>«2),and (18) 
V2to 



V 2to 



IP' 



(19) 



The total mass shift is the sum over those due to in- 
dividual spin channels, as in equation (8). The sum 
over the first term above is trivial on account of equa- 
tion (10). The sum over the second can be done using 
the coefhcients in Table II, and the results are also 
shown there; remarkably the dependence on S and J 
of the sum is proportional to the matrix element of 
the spin-orbit operator. 



Cs'Lf=(l' 



L-s; 



SLJ 



(20) 



with = +1/2 and = ^1/4- Thus the mass shift 
can be written 

AMsLj = J2 (^^'^' + AsLjP' + S{L ■ S)sLjiiY') 

(21) 

The correction to the hypcrfinc splitting due to loops 
follows immediately; everything cancels except a term 
proportional to As, 

- i (AM3p„ 4- 3AM3p^ + SAMapJ + AMip^ 

= -As^P'. (22) 

Thus to this order, in the nonrelativistic limit (Ag = 
0) the result of zero hyperfine splitting survives loop 
corrections. Away from this limit the hyperfine split- 
ting, which is in any case inherently small, is actually 
reduced with respect to its quenched value. In most of 
the quoted examples in Table I, there remains a small 
hyperfine splitting after loop corrections; this is due 
to quadratic corrections to the expansion (16). The 
smallness of -^5^} explains why the mechanism works 
even better for bb than cc. 

In equation (21), the first term is dominant and 
sets the scale of the mass shifts. The smallness of the 
latter two terms is a generic feature of any reasonable 
model, since P' is small, and AM',e << 2to so that 
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Y is also small. In the loop-induced spin-splitting 
between any two states, 

AM S'L'J' - AMsLj, (23) 

the large term cancels, which underlines the observa- 
tion that relative mass shifts are much smaller than 
overall mass shifts [7, 11]. 

Thus there is a hierarchy of scales in the problem. 
While overall mass shifts due to loops can be large, 
the induced spin splittings are small, and the induced 
hyperfine splittings smaller still. The quark model 
prediction of small hyperfine splitting is thus robust 
against corrections due to unquenching. 

If one repeats the entire exercise for D-wave and 
higher L mesons one finds that the same mechanism 
applies. Thus the prediction [15] for the mass of the 
bottomonium in terms of the recently discovered 
[16] ^Di_2.3 ought to be reliable. 

The mechanism depends critically on the coeffi- 
cients Cg^Lj. These are common to all models in 
which quark spin and spatial degrees of freedom fac- 
torise, the coupling is driven by the creation of a 
spin-one pair, and the spin degrees of freedom are 
conserved. The observed small hyperfine splittings 
can thus be interpreted as support for this dynamical 
picture, which has already some support from lattice 
QCD calculations of strong decay [17]. 

Table I suggests two other general properties of 
loop-induced mass shifts. Firstly, the induced hyper- 
fine splitting is always negative. If the physical hyper- 
fine splitting is positive, as is favoured by the bulk of 
experimental and lattice data, then in the absence of 
some other effect the only possibility is that the po- 
tential model splitting is positive. This may help 
to distinguish among different models, which disagree 
on the sign of Ag [18]. 

Secondly, the pattern of mass shifts is always the 
same, 

AAf3p^ > AMip^ > AAhp^ > AMsp^, (24) 

and implies the effect of unquenching is to bring the 
masses closer together with respect to their potential 
model values. (It turns out that this is not unique to 
P-wave levels.) The calculations of Eichten et al. [19] 
exhibit a different pattern of splittings, due to mass 
shifts from mixing between different canonical config- 
urations. Although in their case the induced hyperfine 
splitting remains small, this is not necessarily true of 
all models with configuration mixing. 

Note also that the calculations of Pennington and 
Wilson [20] exhibit a large induced hyperfine splitting; 
their mass shifts are not protected by the mechanism 
outlined here because they do not sum over all spin 
combinations in the loop. The same applies to the 
simplified approach of Shmatikov [21]. 



The expression (21) implies that the effect of un- 
quenching can be absorbed into a renormalisation of 
M, As, At and Ao, and this may have interesting 
phenomenological consequences. For example. At and 
Ao are directly related to the relative contribution 
of vector and scalar parts of the interquark potential 
[21, 22]. 

The factorisation of the angular momentum depen- 
dence of the mass shift by means of the expansion (15) 
is similar in spirit to the approach of Tornqvist and 
Zenczykowski [23], although the expansion parame- 
ters differ. The advantage of the choice in this paper 
is that the first two terms in the expansion involve the 
average mass shift and continuum probability, both of 
which are physically meaningful quantities. The latter 
influences various hadron transition properties, such 
as radiative and pion decay widths. It may be pos- 
sible, by means of the formalism presented here, to 
relate these properties the pattern of spin splittings, 
without making reference to any particular model. 

The author thanks F. Close for useful discussions 
and C. Thomas for pointing out some lattice QCD 
references. 
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